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' Abstract. In the past years a wealth of observations allowed to un- 

ravel the structural properties of the Dark Matter Halos around spirals. 
First, their rotation curves follow an Universal profile (URC) that can 
be described in terms of an exponential thin stellar disk and a dark halo 
^ . with a constant density core, whose relative importance increases with 

• galaxy luminosity. Careful studies of individual objects, from dwarfs 

I to giants, reveal that dark halos have a core, whose size ro correlates 

with the central density po- These properties are in serious discrepancy 
with the cuspy density distribution predicted by N-body simulations 
in coUisionless ACDM Cosmology. 
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Q , Rotation curves (RC's) of disk galaxies are the best probe for Dark Matter (DM) on 
$-H galactic scale since its discovery. However, only recently we discovered some crucial 
c/2 , aspects of their mass distribution by a large number of high-quality RC's and by 
' improvements in the techniques of the RC mass-modeling. The DM distribution 
is usually assumed according to one of following different approaches. 

An empirically one (Persic, Salucci & Stel, 1996, PSS) adopts the simplest halo 
velocity profile that (in combination with the stellar disk) reproduces the Universal 
Rotation Curve of Spirals (out to 3-4 Rd, Ropt = 3.2 disk scale-lenghts Rd)- 

VluRci^) = Volt (1 -/?)(! + , 2^ 2, (1-1) 

where x = R/Ropt, a is the halo velocity core radius, in units of Ropt and /3 = 
(Vd/Vopi)^ is the fractional contribution to the circular velocity of the stellar disk at 
Ropt- This URC-based profile has the advantage of simplicity and, with suitable 
choices for (3 and a, it can represent a variety of mass models, including the 
NFW one. Of course, it cannot be extrapolated beyond the region traced by the 
kinematics. 
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The cosmological approach relies on high-resolution N-body simulations, ac- 
cording to which Cold Dark Matter (CDM) halos achieve the equilibrium den- 
sity profile (Navarro, Frank % White, 1996, NFW): pnfw(-R) = (R/r,){i+R/rs)^ 
where and ps are the characteristic inner radius and density, usually expressed 
in terms of virial mass and radius Tvir, M^r- By setting c = r^ir/rg, since 
from simulations: c ~ 21(M„ir/(10"MQ))-o i3 (Bullock et al. 2001), we have: 
Vnfw = V^FW {R, Myir,c{Myir)), and in detail: 

^NFw(^) — K^ir^^^-j ^ ^ (1-2) 

where x = R/vs and A{x) = ln(l + x) — x/{l + x) 

The third approach adopts the Burkert profile in order to account for the 
observational evidence at inner radii and to converge to the NFW profile at outer 
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Fig. 2. Central halo density po (in g/cm'^) vs. disk mass (in solar units); bottom) vs. 
core radii (in kpc). Details in Salucci & Burkert, 2001 

radii (Salucci & Burkert, 2000): pb{R) = ^?/i°2i i\ with po and being the 
central DM density and the scale radius. Then (Mq — 4 l.Gpo ^o): 

V^J(i?) = G^{ln(l -Fi?/ro) - arctan(i?/ro) + 0.5 ln[l + [R/rof]} (1.3) 

It is important to stress the following points: a) the mass in spirals is dis- 
tributed according to the Inner Baryon Dominance (IBD) regime (PSS, Salucci 
and Persic, 1999): there exists a "transition radius" Ribd — 2i?c; M^^, with 
Md the disk mass in units of 10"Mq. For R < Rjbd, (Salucci et al 2000), the 
luminous matter, in an exponential thin disk of length-scale Ru, accounts for the 
whole gravitating mass. For R > Ribd, instead, an additional dark component, 
distributed unlike the stellar disk, emerges and it rapidly becomes dominant, b) 
the HI and the bulge contribution to V{R) do not play a relevant role in the 
subject of this paper. 

2 Dark Halos Properties from the Universal Rotation Curve 

PSS have derived, from ~ 20000 velocity measurements, relative to ^ 900 rotation 
curves, Vsynij^^] Mj), the synthetic circular velocity of spirals, binned in intervals 
of I- magnitudes (see Fig. 1). At a fixed luminosity and normalized radius, the 
spiral RC's show a cosmic variance, with respect to Vgyn, that is much smaller 
than their radial variations and than their luminosity dependence. As result. 
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spirals sweep a very narrow locus in the RC- profile/amplitude/luminosity space. 
The whole set of synthetic RC's define the Universal Rotation Curve (URC), that 
we analytically represent as the sum of two terms: a) the exponential thin disk: 

V,f{x) = 1.28 P Kp* a:^ {IoKq - /ii^i)|i.6x (2.1) 

and b) the spherical halo given by eq (1.1). Then: V^j^q{x) — V^jjj^(j{x,(5,a) + 
V^jjf^Q{x,j3), with a and (3 as free parameters, reproduces the synthetic curves 
Vsyn{R) up to their rms (i.e. within 2%) when (3 = f3{\ogVopt) and a — a(/3) as 
given in PSS and in Salucci and Burkert, 2000. Notice that also (bulge-free) indi- 
vidual high-quality RC are generally well fit by eq. (l-l)- Inside i?opt smaller ob- 
jects have larger dark-to-steUar mass ratio: M*/Mvir =i 0.2 (M*/2 x lO^^M©)" '^^. 
Scaling relationships relate the halo and disk mass parameters implying that the 
densest halos harbor the least massive disks and the most inefficient ones to trans- 
form the original HI content in present day stars. Notice that the velocity core 
radius a suggests, but docs not prove, a flat core in the DM density distribution, 
that can be unambiguously revealed only by proper individual RC's. 

3 The DM halo density 

The mass structure of spirals is well probed in objects with both HI and H^ high- 
quality high- resolution RC. The existence of the URC and the Universal properties 
of cosmological halos allows us to concentrate only in a reasonable number of 
test cases. The first " absolutely safe" determination of DM halos density profiles 
in Spirals was obtained in Gentile ct al 2004, where 5 spirals with Hi and Ha 
RC's were studied by means of "state of the art" observational, data analysis 
and modeling techniques. In each object HI and Ha rotation curves agree very 
well (where they coexist) and the combined Ha + Hi RC is smooth, symmetric 
and extended out to 6-8 disk length-scales. The mass distribution {\.e.V{R)) is 
modeled as the sum of three components: two stellar /gaseous disks and a spherical 
dark halo: V^^^^^i = Vj^^/, + V| + V^^^, with the halo contribution represented by 
the Burkert profile given by eq (1.3). Light traces the stellar mass via a radially 
constant mass-to-light ratio. The gas contribution Vgas{R) is obtained from HI 
surface brightness and the distance of the object. For each galaxy, we determine the 
values of the structural parameters /?, ro, po by means of a x^-minimization of the 
velocity model: V^^^^i{R; po, P, ro) = V^iR; (3) + V^{R; po, f3, ro) + V^l^iR) to the 
(measured) circular velocity V{R), subject to the constraints: Vmodeiiat Ropt) = 
V{Ropt) and \dV/dR - dVmodei/dR\ < 0.l\dV/dR\ (see Gentile et al. 2004 for 
details). This halo profile (-|- the exp disk) fits the rotation curves extremely well, 
with no systematic deviation. None of the ^ 100 data points of the five RC's 
is discrepant at the 3 a level {a is the observational error). The stellar I-band 
mass-to- light ratios lie between 0.5 and 1.8 and are consistent with population 

synthesis models. The presence of cores is clear: tq = (0.7 2.3) x Rgp^ and 

Po = (0.4 3) X 10^^%cm^'^. It is worth noticing the existence of the relationship 

among the halo structural parameters po = 5 x 10~^^r^^^^ ea;p[— (ro/27)^] g/cm^ 
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Fig. 3. Best-fit models (solid line). Also shown the DM (long- dashed), stellar (dotted) 
and HI (short-dashed) contributions. 

(ro in kpc, Salucci and Burkert, 2001), crucial for its implications on the nature 
of dark matter . 



4 High-quality RC's and NFW halos 

The same objects, when modeled by means of the stellar/gas disks + NFW halo 
components fail to reproduce the shape of the observed rotation curve. Moreover, 
they show systematic discrepancy in the predicted velocities: these, in the central 
parts, are too high. In detail: 10% of the measurements cannot be matched in 
any way, the difference between them and the predictions exceeding Scr (i.e. three 
times the observational error). An other 10% of data suffer of a poor match, the 
offset is at the level of 2-3 a. This is a huge model-data discrepancy: since 1980, 
every mass model without a central cusp has predicted circular velocities that 
resulted within, at worst, 2 a the observed ones. 

Let us notice that, also if we leave c as a free parameter, there is no appreciable 
improvement in the model fits, and that, at galactic scales, the actual value of the 
density inner slope is around -1.3, (Navarro et al 2004) making things even more 
difficult for standard ACDM. A further shortcoming is that the resulting NFW 
disk mass-to-light ratio turn out to be unaccoptably lower than the values we 
estimate from the galaxy color (Gentile et al, 2004). 

Finally, we draw attention on the further evidence provided by de Blok and 
Bosma, 2002, Simon et al 2003, Wcldrakc et al, 2003, Bolatto et al 2002 about 
the serious theory vs observations discrepancy in the DM density distribution, we 
show in Fig 3 a test case. 
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Fig. 4. The density of the dark halo of 116-G12,and the CDM predictions right 

5 Conclusions 

The mass distribution in spirals underlies tight relationships among their struc- 
tural quantities, likely, as the results of strong feedbacks occurred during early 
stages of galaxy formation. DM halos around galaxies have an inner constant- 
density region, whose size exceeds the stellar disk length-scale and emerge as an 
one-parameter family. The order parameter (either the central density or the core 
radius) correlates with the stellar mass. There is no evidence that the density 
profile converges, at large radii, to a p ~ (or steeper) profile. The DM dis- 
tribution is determined by physical parameters, the central core density and the 
core radius, that have no counterpart in the gravitational instability/hierarchical 
clustering picture. 

Solutions for the existence of a region of "constant" density include a) DM 
"interacted" with baryons. Original "cuspy" halos have been smoothed out b) DM 
has a different power spectrum/perturbations evolution than the current Standard 
Picture c) Dark Matter is a "field", that mimics the effects of a cored halo of 
particles, d) the actual dynamical evolution of DM halos, including their baryonic 
content, is more complex than that presently emerging in simulations. 

References 

J.S. Bullock, et al MNRAS, 321, 559 (2001) 
W.J.G. de Blok, A. Bosma, AA, 385, 816 (2002) 
G. Gentile et al. MN 351, 903 (2004) 

J.F. Navarro, C.S. Frenk, S.D.M. White: ApJ, 462, 563 (1996) 

J.F. Navarro et al : MN, 349, 1039 (2004) 

M. Persic, P. Salucci, F. Stel: MNRAS, 281, 27P (1996) 

P. Salucci, M. Persic: AA, 351, 442 (1999) 

P. Salucci, A. Burkert: ApJ, 537L, 9S (2000) 

A. D. Bolatto et al ApJ, 565, 238 (2002) 

J. D. Simon et al ApJ, 596, 957 (2003) 

D. Weldrake et al., MN 340, 12 (2003) 



Dark Matter in Spiral Galaxies 

G. Mamon, F. Combes, C. Deffayet, B. Fort (eds) 

EAS Publications Series, Vol. ?, 2005 



THE DISTRIBUTION OF DARK MATTER IN SPIRALS 

Paolo Salucci^ 

Abstract. In the past years a wealth of observations allowed to un- 
ravel the structural properties of the Dark Matter Halos around spirals. 
First, their rotation curves follow an Universal profile (URC) that can 
be described in terms of an exponential thin stellar disk and a dark halo 
with a constant density core, whose relative importance increases with 
galaxy luminosity. Careful studies of individual objects, from dwarfs 
to giants, reveal that dark halos have a core, whose size ro correlates 
with the central density po- These properties are in serious discrepancy 
with the cuspy density distribution predicted by N-body simulations 
in collisionless ACDM Cosmology. 

1 Introduction 

Rotation curves (RC's) of disk galaxies are the best probe for Dark Matter (DM) on 
galactic scale since its discovery. However, only recently we discovered some crucial 
aspects of their mass distribution by a large number of high-quality RC's and by 
improvements in the techniques of the RC mass-modeling. The DM distribution 
is usually assumed according to one of following different approaches. 

An empirically one (Persic, Salucci & Stel, 1996, PSS) adopts the simplest halo 
velocity profile that (in combination with the stellar disk) reproduces the Universal 
Rotation Curve of Spirals (out to 3-4 Rd, Ropt = 3.2 disk scale-lenghts Rd)'- 

V^uRci^) = Kit (1 -/?)(! + a') (^2^ (1-1) 

where x = R/Ropt, a is the halo velocity core radius, in units of Ropt and /3 = 
{Vd/Vopt)'^ is the fractional contribution to the circular velocity of the stellar disk at 
Ropt- This URC-based profile has the advantage of simplicity and, with suitable 
choices for /? and a, it can represent a variety of mass models, including the 
NEW one. Of course, it cannot be extrapolated beyond the region traced by the 
kinematics. 
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Fig. 1. Synthetic rotation curves (circles) and tlie Universal Rotation Curve (solid line). 
The dark/luminous contributions are indicated with a dashed/dotted line. 



The cosmological approach relics on high-resolution N body sinnilations, ac- 
cording to which Cold Dark Matter (CDM) halos achieve the equilibrium den- 
sity profile (Navarro, Prenk % White, 1996, NFW): /9nfw(-R) = (fl/r,)(i+fl/r,)a 
where r,, and pa are the characteristic inner radius and density, usually expressed 
in terms of virial mass and radius fvir, M^ir- By setting c = r^ir/rs, since 
from simulations: c ~ 21{Myir/{10^^MQ))-°-'^^ (Bullock et al. 2001), we have: 
Vnfw = V^FW {R, Myir,c{Myir)), and in detail: 

^NFw(-^) = K^ir^^^^ ^ (1-2) 

where x = R/rg and A{x) = ln(l + x) — x/{l + x) 

The third approach adopts the Burkert profile in order to account for the 
observational evidence at inner radii and to converge to the NFW profile at outer 
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radii (Salucci & Burkert, 2000): pb{R) — (j;^rp°(jj'2_)_^ii) with po and being the 
central DM density and the scale radius. Then (A/q = 4 1.6po ''o): 

V^{R) = G^{ln(l + i?/ro) - arctan(i?/ro) + 0.5 ln[l + [R/rof]} (1.3) 

It is important to stress the following points: a) the mass in spirals is dis- 
tributed according to the Inner Baryon Dominance (IBD) regime (PSS, Salucci 
and Persic, 1999): there exists a "transition radius" Ribd — '^Rd ^y^, with 
Md the disk mass in units of 1O^^M0. For R < Ribd, (Salucci et al 2000), the 
luminous matter, in an exponential thin disk of length-scale Rb, accounts for the 
whole gravitating mass. For R > Ribd, instead, an additional dark component, 
distributed unlike the stellar disk, emerges and it rapidly becomes dominant. 6) 
the HI and the bulge contribution to V{R) do not play a relevant role in the 
subject of this paper. 

2 Dark Halos Properties from the Universal Rotation Curve 

PSS have derived, from ^ 20000 velocity measurements, relative to ^ 900 rotation 
curves, Vsyn ( ^ i Mi), the synthetic circular velocity of spirals, binned in intervals 
of I- magnitudes (see Fig. 1). At a fixed luminosity and normalized radius, the 
spiral RC's show a cosmic variance, with respect to Vsyn, that is much smaller 
than their radial variations and than their luminosity dependence. As result, 
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spirals sweep a very narrow locus in the RC- profile/amplitude/luminosity space. 
The whole set of synthetic RC's define the Universal Rotation Curve (URC), that 
we analytically represent as the sum of two terms: a) the exponential thin disk: 

V^{x) = 1.28 /? {loKo - /ii^i)|i.6x (2.1) 

and b) the spherical halo given by eq (1.1). Then: V^i^q{x) = jjj^(-,{x, (5,a) + 
V^ifj^Q{x,P), with a and (3 as free parameters, reproduces the synthetic curves 
Vayn{R) up to their rms (i.e. within 2%) when /3 = /3(log Vopt) and a = a{(3) as 
given in PSS and in Salucci and Burkert, 2000. Notice that also (bulge-free) indi- 
vidual high-quality RC are generally well fit by eq. (1.1). Inside i?opt smaller ob- 
jects have larger dark-to-stellar mass ratio: M*/Mvir ^ 0.2 (M*/2 x lO^Mo)^-'^^. 
Scaling relationships relate the halo and disk mass parameters implying that the 
densest halos harbor the least massive disks and the most inefficient ones to trans- 
form the original HI content in present day stars. Notice that the velocity core 
radius a suggests, but does not prove, a flat core in the DM density distribution, 
that can be unambiguously revealed only by proper individual RC's. 

3 The DM halo density 

The mass structure of spirals is well probed in objects with both HI and H^ high- 
quality high-resolution RC. The existence of the URC and the Universal properties 

of cosmological halos allows us to concentrate only in a reasonable number of 
test cases. The first "absolutely safe" determination of DM halos density profiles 
in Spirals was obtained in Gentile et al 2004, where 5 spirals with Hi and Ha 
RC's were studied by means of "state of the art" observational, data analysis 
and modeling techniques. In each object HI and Ha rotation curves agree very 
well (where they coexist) and the combined Ha -|- Hi RC is smooth, symmetric 
and extended out to 6-8 disk length-scales. The mass distribution {i.e.V{R)) is 
modeled as the sum of three components: two stellar/gaseous disks and a spherical 
dark halo: V^^^^i = Vj^j, + -|- V^^^, with the halo contribution represented by 
the Burkert profile given by eq (1.3). Light traces the stellar mass via a radially 
constant mass-to-light ratio. The gas contribution Vgas{R) is obtained from HI 
surface brightness and the distance of the object. For each galaxy, we determine the 
values of the structural parameters /?, rg, po by means of a -minimization of the 
velocity model: Vl^^^i{R- po, l3, ro) = V^(R- /?) + V|(i?; po, ^, tq) + V^^^{R) to the 
(mcasmcd) circular velocity V{R), subject to the constraints: Vmodei{at Ropt) = 
V{Ropl) and \dV/dR - dVmodei/dR\ < 0.1\dV/dR\ (sec Gentile et al. 2004 for 
details). This halo profile (+ the exp disk) fits the rotation curves extremely well, 
with no systematic deviation. None of the ~ 100 data points of the five RC's 
is discrepant at the 3 a Icvoi {a is the observational error). The stellar I-band 
mass-to-light ratios lie between 0.5 and 1.8 and are consistent with population 

synthesis models. The presence of cores is clear: ro = (0.7 2.3) x Ropt and 

po = (0.4 3) X lO^^'^gcmr^ . It is worth noticing the existence of the relationship 

among the halo structural parameters po = 5 x 10~-^''rg '^^'^ ea;p[— (ro/27)^] g/cm^ 
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Fig. 3. Best-fit models (solid line). Also shown the DM (long- dashed), stellar (dotted) 
and HI (short-dashed) contributions. 

(ro in kpc, Salucci and Burkert, 2001), crucial for its implications on the nature 
of dark matter . 



4 High-quality PC's and NFW halos 

The same objects, when modeled by means of the stellar/gas disks + NFW halo 
components fail to reproduce the shape of the observed rotation curve. Moreover, 
they show systematic discrepancy in the predicted velocities: these, in the central 
parts, are too high. In detail: 10% of the measurements cannot be matched in 
any way, the diffcrc;nc;c; between them and the predictions exceeding 3a (i.e. three 
times the observational error). An other 10% of data suffer of a poor match, the 
offset is at the level of 2-3 a. This is a huge model-data discrepancy: since 1980, 
every mass model without a central cusp has predicted circular velocities that 
resulted within, at worst, 2 a the observed ones. 

Let us notice that, also if we leave c as a free parameter, there is no appreciable 
improvement in the model fits, and that, at galactic scales, the actual value of the 
density inner slope is around -1.3, (Navarro et al 2004) making things even more 
difficult for standard ACDM. A further shortcoming is that the resulting NFW 
disk mass-to-light ratio turn out to be unacceptably lower than the values we 
estimate from the galaxy color (Gentile et al, 2004). 

Finally, we draw attention on the further evidence provided by de Blok and 
Bosma, 2002, Simon et al 2003, Weldrake et al, 2003, Bolatto et al 2002 about 
the serious theory vs observations discrepancy in the DM density distribution, we 
show in Fig 3 a test case. 
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Fig. 4. The density of the dark halo of 116-G12,and the CDM predictions right 

5 Conclusions 

The mass distribution in spirals underlies tight relationships among their struc- 
tural quantities, likely, as the results of strong feedbacks occurred during early 
stages of galaxy formation. DM halos around galaxies have an inner constant- 
density region, whose size exceeds the stellar disk length-scale and emerge as an 
one-parameter family. The order parameter (either the central density or the core 
radius) correlates with the stellar mass. There is no evidence that the density 
profile converges, at large radii, to a p ^ (or steeper) profile. The DM dis- 
tribution is determined by physical parameters, the central core density and the 
core radius, that have no counterpart in the gravitational instability/hierarchical 
clustering picture. 

Solutions for the existence of a region of "constant" density include a) DM 
"interacted" with baryons. Original "cuspy" halos have been smoothed out b) DM 
has a different power spectrum/perturbations evolution than the current Standard 
Picture c) Dark Matter is a "field", that mimics the effects of a cored halo of 
particles, d) the actual dynamical evolution of DM halos, including their baryonic 
content, is more complex than that presently emerging in simulations. 
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